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We present a broad class of states which are diagonal in the basis of A^-qubit GHZ states such that 
non-positivity under the partial transpose operation is necessary and sufficient for the presence of 
entanglement. This class includes many naturally arising instances such as dephased or depolarised 
GHZ states. Furthermore, our proof directly leads to an entanglement witness which saturates 
this bound. The witness is applied to thermal GHZ states to prove that the entanglement can be 
extremely robust to system imperfections. 



Introduction: Multipartite entanglement is still a phe- 
nomenon that is poorly understood and categorised. For 
some types of entangled state, such as 7V-qubit GHZ 
states, it requires very little noise (loss of a single qubit) 
to entirely destroy the entanglement, whereas others, 
such as the multipartite states based on error correct- 
ing codes, seem much more robust. If we are ever to use 
entanglement as a resource in information processing pro- 
tocols, it is essential to understand when entanglement is 
present in a system and how to detect it. For instance, 
the ability of entanglement to persist at high tempera- 
tures could vastly reduce experimental requirements, and 
has a direct bearing on the possibility of constructing 
quantum memories [1 . 

Attempts to characterise and detect multipartite en- 
tanglement are certainly not new. See, for instance, |2j- 
[lOj and, in particular, the wide-ranging review of [9] and 
the further references contained therein. While the ma- 
jority of these strategies are capable of detecting some 
entanglement, the majority, with some notable excep- 
tions including |2H5] , are unable to convey how well these 
characterisations perform. They might be capable of de- 
tecting some entanglement, but how much is missed? 

It seems a reasonable starting point for these studies is 
to consider classes of states, such as those that are diago- 
nal in the GHZ basis, which arise frequently in quantum 
information, and are hence likely to be of most inter- 
est. Some partial categorisations are already known. For 
instance, thermal GHZ states can be distilled up to a fi- 
nite temperature and this temperature is tight i.e. above 
that temperature, entanglement can't be distilled [4, 5 . 
There is still entanglement present in these models above 
the distillation threshold [6] , which is consequently bound 
entanglement. Diir and Cirac [3] also considered a sub- 
set of GHZ-diagonal states and showed that a necessary 
and sufficient condition for distillation of this subset is 
that the state should be non-positive with respect to the 
partial transpose (NPT) operation across every possible 
bipartition, whereas a necessary and sufficient condition 
for full separability is that the state should be positive 
with respect to the partial transpose (PPT) across all 
possible bipartitions. It does not seem altogether sur- 
prising that these two conditions do not, in general, co- 
incide, and hence it is interesting to understand, in as 



wide a context as possible, how persistent entanglement 
can be, and how to detect it. 

In this paper, we give a broad class of GHZ-diagonal 
states for which the PPT condition is necessary and suf- 
ficient for the state to be fully separable. This class in- 
cludes the special cases of the thermal state studied in 
[IHS] and the GHZ-diagonal states in [3]. Our formal- 
ism instantly yields an entanglement witness for opti- 
mally detecting that entanglement, as well as providing 
the foundation for a future, more wide-ranging, study of 
identical concepts in thermal graph states [11 . We de- 
scribe the special cases of the thermal GHZ state (which 
is also a dephased GHZ state) and a depolarised GHZ 
state. 

GHZ-diagonal states: Consider a system of N qubits, 
and define the stabilizers 

\ ZiXn n>2 

These terms mutually commute. is the Pauli X ma- 
trix applied to qubit n, and Z^ denotes the application 
of Z rotations to all qubits for which the A/'-bit string x 
is 1, i.e. by denoting the n^^ bit of x G {0, 1}^ as X Yi 5 

N 
n=l 

Kx is similarly defined in terms of the K^. The +1 eigen- 
state of each of these stabilizers is 

|^) = |0)|+)«(^-i) + |l)|-)«(^-^), 

and = Zx\ip) for x G {0, 1}^ is an eigenstate of all 
products of stabilizers Ky with eigenvalue (—1)^"^. There 
is a straightforward local equivalence to any other GHZ 
basis that one might choose, which does not affect the 
entanglement structure that we are investigating. How- 
ever, this particular formulation will allow an immediate 
translation of many of our results to the more general 
case of graph states. Any state 
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is diagonal in this basis, and it is this class of states which 
we consider. In order for p to be a valid state, we require 
5o = 1 and min^^jo,!}^ ^y{~^Y'^ ^ (i.e. the eigen- 
vector with minimum eigenvalue is non- negative). 

We are now interested in evaluating the partial trans- 
pose criterion on this state in a first step to determine 
when there is entanglement present in the state. Start- 
ing from Eqn. ([T]), we introduce a bipartition z G {0, 1}^ 
i.e. all the vertices with = are on one side of the 
partition and those with = 1 are on the other side, 
and take the partial transpose on the Zn = I side. With- 
out loss of generality, qubit 1 can be placed on the 
side of the bipartition and hence z is restricted to being 
just over the N — 1 other qubits. Recall that under the 
partial transpose, the Pauli operators alter by Z^, 

Xn ^ Xn but Yn ^ {-ly-'Yn. ThuS, 
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Observe that products of stabilizers remain as products 
of stabilizers and, as a result, the eigenvectors of p^^ are 
just IV^a^), with eigenvalues fx^z{s)/2^: 



fx, 



ye{0,l}' 



E^ 



If there exists a choice of x and z such that fx,z{s) < 0, 
the state is certainly entangled due to being NPT across 
the bipartition z. 

Entanglement Witnesses Saturating PPT: Using this 
formalism, it's straightforward to find an entanglement 
witness that will saturate the PPT threshold for any state 
which is diagonal in the GHZ state basis. To do this, we 
measure the observables 

2yG{0,l}^ 
aG{0,l} 



a{x®z) I 



For any arbitrary density matrix p with GHZ stabilizer 
expectation values 5, 

Tt{Wx,,p)= {-irn-^y'^-''-'-'Sy = fx^,{s). 
ye{o,i}^ 

Hence, for GHZ-diagonal states, this gives the eigenval- 
ues of the partial transpose of the state about bipartition 
2:, and finding Ti (Wx^zp) < for any x or z proves it's en- 
tangled. This is a genuine entanglement witness in that, 
for any state p = Y^^.y l^^.v \^x) {'^y\, which may not be 
diagonal in the graph state basis, finding one of the ob- 
servables to be negative witnesses the fact that it's entan- 
gled. To prove this, note that any p can be converted, via 
local probabilistic operations, into a graph diagonal state 
Pd = Yxl^x.x \'^x) (V^xl with the same diagonal elements 
[12], and hence the same values of s. So, if p is fully 



separable, it will have the same value of Tr{W p) as pd-, 
which we know will be positive since the local conversion 
to a diagonal state cannot introduce entanglement. 

Separability: We are now in a position where we can 
determine whether a GHZ-diagonal state is NPT with 
respect to some bipartition, and have an observable that 
can witness the entanglement. We now move to studying 
the converse, when the state is certainly not entangled. 
Again, the stabilizer formalism is immensely helpful. We 
will say that Kx and Ky have a compatible basis if at ev- 
ery site n when Kx is a Pauli matrix a, then Ky is either 
<j or 1 at that site and vice versa. We are interested in 
such cases because each product of stabilizers Kx is just 
a tensor product of Pauli operators, and hence its eigen- 
vectors are product states. Two terms Kx and Ky have 
a simultaneous product state decomposition if they have 
a compatible basis. So, in order to give a fully separa- 
ble decomposition of p, we can simply group together all 
terms that have a compatible basis, and find the small- 
est eigenvalue. This grouping of terms has to have some 
component of the 1 added such that the minimum eigen- 
value is 0. If we do this, then that grouping of terms is 
a separable state, with a decomposition specified by the 
common product basis. We are finally left with a condi- 
tion that the excess weight of 1 terms should be positive. 
In the case of a GHZ-diagonal state, the terms Ky for 
y G {0, 1}^ with ^1 = 1 do not have any compatible 
terms, whereas all terms Ky with ^1 = are mutually 
compatible. We hence change notation slightly to Kiy 
and K^y respectively for y G {0, The decomposi- 
tion therefore takes the form 



E 

ye{o,i}'' 



+1 (.,,-^„_.E«o.(-ir--ENi.i)- 



Thus, provided 



min ^o,(-ir^- Y l^i.ll>0' 

2yG{0,l}^-i ^e{0,l}^-i 



we have a separable decomposition of p. Compare this 
to fxox,x®~z{^) where (-1)^° = -sgn(5ioo...o), 

Y ^o,(-l)'-^-sgn(.ioo...o) Y ^i.(-l)'"'- 

yG{0,l}^-i 2yG{0,l}^-i 

The two are equal if x corresponds to the minimal choice 
in the separable state decomposition, and there exists a 
z e {0, 1}^-^ such that 



sgn(5ioo...o)5iy(-l)"^ > G {0, 1} 



AT-l 



(2) 



(If S100...0 = then xq remains a free parameter.) If this 
simple condition is satisfied, then PPT exactly detects 
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FIG. 1: Comparison of the PPT critical temperature for the 
thermal GHZ state, a lower bound for the model when per- 
turbed by a uniform magnetic field of strength 5 / /\ — 0.3, and 
performance comparison of a previous entanglement witness 
[6]. Choice ofA = A:s = ^= l ensures unitless quantities. 

the transition between the existence of bipartite entan- 
glement and full separability of the state which, in turn, 
makes our entanglement witnesses optimal. 

Eqn. ([2| gives a sufficient condition for the coincidence 
of thresholds for full separability and PPT. If not satis- 
fied, is there really a separation between the PPT thresh- 
old and the best known separable state? For a 3-qubit 
GHZ-diagonal state, Eqn. Q is fulfilled provided 



n 

2yG{0,l}2 



> 0. 



We conclude that roughly half of the parameter space is 
covered by Eqn. ([2| in this case. One example outside 
this regime is the state 



]\{l^Kn)-2KiK^^al 



8(1 + a) 



Provided a > 2, p is a valid state, but it is also PPT 
with respect to all possible bipartitions. Our previous 
construction of a separable state is valid for a > 4. This 
can be improved to a > 2^/2 by rewriting the sum Ki + 
K1K2 - KiKs + K1K2KS as 

1 

i ^(x + (-irr)i(^ + {-irY)2{z - (-i)"y)3. 

n=0 

Upon implementing the semi-definite programming tech- 
niques of [13 , we witnessed entanglement in the region 
a < 2.828. Hence, the separable decomposition is not 
universally optimal, but neither is the PPT condition. 

Thermal States and Perturbations: Many noise mod- 
els satisfy Eqn. ([2]), including the sub-class considered in 
[3] (all Sly equal and positive). We will now discuss two 
special cases. The first is local dephasing noise on each 
qubit with probability p, and the second will be depo- 
larising noise. Dephasing noise also corresponds to the 



thermal state of the Hamiltonian 



N 



which has Sy = Yl'^^-^tSinh{l3An/2)y 
Hamming weight of /3 is the inverse temperature, and 



where Wy is the 



are energy terms and p = e" )~ . Since Sy 
z = 00 ... 0. One can also check that x is 11 ... 1 
consequently derive a simple threshold condition. 



>0, 
and 



tanh(/3Ai/2) = e-^^^=2A.^ 

Even though distillable entanglement only persists to a 
finite temperature 01 15], bound entanglement, which is 
all bipartite and detected by the partial transpose, per- 
sists to a temperature that increases with N (see Fig.jl]). 
Previous entanglement witnesses, even those specifically 
designed to detect GHZ state entanglement [6], were of- 
ten far from optimal. As observed numerically for other 
graph states in [14 , the entanglement is very robust to 
perturbations in the A^. We can also add a local mag- 
netic field term 



N 



N 



n=l 



n=l 



Since the terms K'^ = AnKn + SnZn mutually commute, 
1 



For simplicity of notation, take all the A^ equal, and all 
the Sn equal, although none of the following depends on 
it. We also set s = tanh(/3VA2 -^S^/2). The thermal 
state can be expanded as 



2^p = 



E 



E 



!:e{0,l}" yG{0,l}^ 



VA2 + 



Wx-\-Wy 



S^^/^^VZ^Ky 



although the summation over y is restricted to cases 
where ^/^ = if = 1- To prove the presence of en- 
tanglement, we can use the entanglement witness Wx,x®z 
from the unperturbed case. 



Tr(V^5,5e.~P) 



rn 1 1 AT V 



-sA 



ye{o,i} 



VA2 + (52 



(-1) 



EN 
n=2 ynZn-1 



since TT{ZxKyKz) = Sx^oo...o^y,z- The critical inverse 
temperature Ps at which the expectation value of this 
state is zero is hence related to the unperturbed /3o by 



VA2 + (52 



tanh(^/35\/A2 + (52) = tanh(^/3oA). 



Furthermore, /3s is an upper bound on the true critical 
/3, i.e. a lower bound on the critical temperature. Fig. [l] 
indicates just how robust this entanglement is. 
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FIG. 2: The increasing tolerance of a GHZ state to dephas- 
ing and depolarising noise with system size. Note that the 
two thresholds should not be directly compared since, in the 
case of depolarising noise, one could argue that an error only 
happens with a probability 3p/4, not p. 



Depolarising noise: GHZ states (in our chosen basis) 
are quite robust against dephasing noise, increasingly so 
with the number of qubits involved. They are also very 
tolerant of different noise rates (A^) on different qubits. 
Another experimentally relevant scenario is when the 
GHZ state is depolarised independently on each qubit. 
This is described by a map (on all qubits n) of 

Snip) = (1 - P)P + I (P + XnpXn + Y^pYn + Z^pZ^) • 

With a little effort, one can prove that siy = (l—p)^ and 
soy = (1 — p)2r^y/2l Since all siy are equal, [3 applies 
for both distillability and separability properties. In the 
regime of < p < 1, s^y > and hence z = 00 ... 0. To 
find X, we note that 

min V i-ir'soy 

^ ^ 2yG{o,i}^-i 

depends only on the weight of x. With this observation 
in place, one can prove that 

^{-ir'^soy = 1(2 + ^(2 

y 



with a minimum of (2 — L^/2JpL^/2J occurring at = 
[N/2\ . Hence, the PPT threshold and full separability 
boundary occurs at 

(2 -p)W2JpW2J =2^-1(1 -p)^. 

The critical p still increases with in a similar fashion 
to the dephasing noise. Fig. |2j 

Conclusions: We have given a sufficient condition, 
which naturally encompasses a vast range of GHZ- 
diagonal states, including those that are experimentally 
relevant, such that the existence of an NPT bipartition 
is necessary and sufficient for the state to be entangled. 
We have also shown that there are examples not covered 
by this condition such that the statement is not true. 
We have described entanglement witnesses that detect 
the existence of an NPT bipartition. These are vastly 
stronger than previous witnesses, and simply correspond 
to measuring the overlap with four different states. This 
witness shows that the entanglement in thermal states 
of GHZ graphs is extremely robust to some classes of 
perturbation. Future work could focus on witnessing en- 
tanglement in those cases where entanglement persists 
outside the PPT regime, beyond the current reliance on 
numerical techniques. Criteria developed in [10 can al- 
ready prove the existence of entanglement beyond the 
PPT threshold. 

We have been careful to express much of this paper 
in very general terms using the stabilizer formalism. As 
such, much of the work extends to graph-diagonal states, 
which are also defined by stabilizers, and which cover 
many of the interesting states in quantum information, 
such as cluster states and error correcting codes. We will 
study this case in more detail in [TT]. 
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